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Abstract.  We  consider  a two-level  quantum  system  (qubit)  which  is  continuously  measured  by 
a detector.  The  conventional  formalism,  which  implies  the  ensemble  averaging,  describes  the 
gradual  decoherence  of  the  qubit  state  due  to  measurement.  However,  in  each  particular  realization 
of  the  measurement  process  we  can  have  the  opposite  effect:  gradual  purification  of  the  qubit 
density  matrix.  This  can  be  described  by  the  recently  developed  Bayesian  formalism  suitable  for 
individual  quantum  systems.  The  purification  effect  may  be  verified  experimentally  using  present- 
day  technology  and  can  be  useful  for  quantum  computing.  In  particular,  the  decoherence  of  a single 
qubit  can  be  suppressed  using  continuous  measurement  and  the  feedback  loop. 


The  significant  progress  in  experimental  techniques  during  recent  years  as  well  as  the 
active  research  on  quantum  computing  have  motivated  renewed  interest  in  the  problems  of 
quantum  measurement,  including  the  long-standing  “philosophical”  questions.  In  contrast 
to  the  usual  case  of  averaging  over  a large  ensemble  of  similar  quantum  systems,  it  is 
becoming  possible  to  study  experimentally  the  evolution  of  an  individual  quantum  system. 
In  this  paper  we  consider  the  continuous  measurement  of  a qubit  (two-level  system)  state 
by  a “weakly  responding”  detector  which  can  be  treated  as  a classical  device  [!].  As 
examples,  we  will  discuss  the  following  solid-state  realizations  of  such  a measurement. 
First,  the  location  of  a single  electron  in  a double-quantum-dot  can  be  measured  by  a 
nearby  quantum  point  contact  (QPC)  in  such  a way  that  the  QPC  barrier  height  and,  hence, 
the  current  through  the  detector  are  sensitive  to  the  measured  electron  position  | , ].  The 
second  possible  setup  is  a single  Cooper  pair  box,  the  charge  state  of  which  is  measured 
by  a capacitively  coupled  single-electron  transistor  (SET)  | ].  Finally,  the  flux  state  of  a 
SQUID  can  be  continuously  measured  by  another  inductively  coupled  SQUID  | ]. 

The  conventional  approach  describes  the  measurement  process  by  the  following  equa- 
tions for  the  qubit  density  matrix  pij  in  the  basis  of  “localized”  states: 

H 

Pit  = -P22  - -2  — Impi2,  (1) 

n 

,e  . H 

Pi2  = i ~P  12  + i — (P li  — P22)  — Tpi2  , (2) 

ri  n 

where  H is  the  mixing  (tunneling)  and  s is  the  energy  asymmetry  of  the  qubit  states,  while  T 
is  the  dephasing  due  to  measurement.  In  the  case  of  the  QPC  as  a detector  T = (AI)2/4Sj 
| , , ] where  A/  = I2  — I\  is  the  difference  between  the  average  detector  currents 
corresponding  to  two  localized  states  of  the  qubit  (we  assume  | A/|  <SC  (/j  + h)/2)  and  5/ 
is  the  low  frequency  spectral  density  of  the  detector  shot  noise. 
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Notice  that  Eqs.  ( I )— (2)  do  not  depend  on  the  detector  output  that  is  a consequence 
of  averaging  over  the  ensemble  of  systems.  The  situation  is  completely  different  in  the 
case  of  an  individual  quantum  system  since  the  system  evolution  should  become  dependent 
(“conditioned”)  on  the  particular  detector  output.  The  theory  of  conditioned  (selective) 
evolution  of  a pure  wavefunction  was  developed  relatively  long  ago,  mainly  for  the  pur- 
poses of  quantum  optics  (see,  e.g.  Ref.  [ ] and  references  therein).  However,  for  solid 
state  structures  the  problem  of  continuous  quantum  measurement  with  an  account  of  the 
measurement  result  has  been  addressed  only  recently  [ ],  with  the  main  emphasis  on  mixed 
quantum  states  and  detector  nonideality. 

In  the  Bayesian  formalism  developed  in  Ref.  | ] the  evolution  of  the  qubit  density  matrix 
p is  described  by  the  equations 

(3) 

A)]  pn  — y P\2,  (4) 

where  I (f ) is  the  particular  detector  output,  Iq  = (l\  + h) /2,  and  y is  the  extra  dephasing 
due  to  the  “pure  environment”, 

y = r — (a/)2/4S7.  (5) 

There  is  no  extra  dephasing,  y = 0,  for  the  measurement  by  a QPC,  which  thus  represents 

an  ideal  detector,  while  the  SET  in  a typical  operation  point  is  a significantly  nonideal 

detector,  y ~ T [ ].  The  SQUID  can  be  an  ideal  detector  only  if  its  total  sensitivity  is 
quantum-limited  [ | (A/2  in  energy  units). 

Equations  (3)-(4)  allow  us  to  calculate  the  evolution  of  the  system  density  matrix  for  a 
given  detector  output  / (f).  In  order  to  analyze  the  behavior  of  I ft),  these  equations  should 
be  supplemented  by  the  formula 

W)-Io=  ^-(P22-Pn)  + Ht),  (6) 

where  the  zero-correlated  (“white”)  random  process  f (f ) has  zero  average  and  the  same 
spectral  density  as  the  detector  current,  Sj  — Sj.  (We  use  the  Stratonovich  formalism  for 
stochastic  differential  equations.)  Notice  that  even  though  the  Bayesian  formalism  is  valid 
only  for  “weakly  responding”  detectors,  | A/|  <$C  /o,  the  dimensionless  coupling  AT ///  is 
arbitrary,  so  the  formalism  can  be  used  in  the  Quantum  Zeno  regime  as  well  as  in  the  case 
of  weakly  perturbed  quantum  (Rabi)  oscillations. 

Figure  1 shows  the  result  of  the  Monte-Carlo  simulation  of  the  continuous  measurement 
by  a slightly  nonideal  detector,  y — 0. 1 E,  in  the  case  when  the  evolution  starts  from  the 
maximally  mixed  state,  p\  \ — P22  — 0.5,  p 12  = 0.  One  can  see  that  p\ 2 gradually  appears 
during  the  measurement,  eventually  leading  to  well-pronounced  quantum  oscillations.  In 
the  case  y = 0 the  density  matrix  becomes  pure  after  a sufficiently  long  time.  This  gradual 
purification  can  be  interpreted  as  being  due  to  the  gradual  acquisition  of  information  about 
the  system.  The  detector  nonideality,  y ^ 0,  causes  decoherence  and  competes  with  the 
purification  due  to  measurement. 

In  contrast  to  QPC,  the  SET  as  a detector  directly  affects  the  qubit  energy  asymmetry  e 
because  of  the  fluctuating  potential  0(f)  of  SET’s  central  island.  Since  there  is  typically  a 
correlation  between  fluctuations  of  I ( f ) and  0(f)  [ ],  Eqs.  (3)-(4)  can  be  further  improved 


Pi  1 = — P22  = — 2 — Impi2 
A 


Pi  1P22  |2(f)  — /0], 


. e . H A/ 

P12  = i 7 P12  + i -T-  (pn  — P22)  + — (Pn  — P22)  1 2(f) 
A A Si 
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Fig.  1.  Gradual  purification  of  the  qubit  density  matrix  p (/ ) in  a course  of  continuous  measurement. 


taking  into  account  the  information  about  the  likely  fluctuations  s(t)  caused  by  the  SET.  It 
is  natural  to  add  into  Eq.  (4)  the  term  ip\2K[I(t)  — (p\\J\  + PiihW  = ipnK^it).  where 
K — {ds/d^Sfj,! /Sjh  and  S<pi  is  the  mutual  low-frequency  spectral  density.  One  can 
easily  check  that  the  addition  of  this  term  corresponds  to  the  partial  recovery  of  coherence, 
so  that  the  dephasing  rate  y should  be  replaced  with  y = y — K2S//A. 

To  observe  the  density  matrix  purification  experimentally,  it  is  necessary  to  record  the 
detector  output  with  sufficiently  wide  bandwidth,  A / T (possibly,  A / ~ 109  Hz),  and 

plug  it  into  Eqs.  (3)-(4).  Calculations  will  show  the  development  of  quantum  oscillations 
with  precisely  known  phase.  Stopping  the  evolution  by  rapidly  raising  the  qubit  barrier 
(H  —*■  0)  when  p\  \ ~ I and  checking  that  the  system  is  really  localized  in  the  first  state, 
it  is  possible  to  verify  the  presented  results. 

Density  matrix  purification  can  be  used  in  quantum  computing  to  suppress  the  gradual 
qubit  decoherence  due  to  interaction  with  the  environment  (to  keep  a qubit  “fresh”).  The 
idea  is  to  use  the  feedback  loop  which  controls  the  qubit  parameter  H (control  of  e is 
also  possible)  in  order  to  decrease  the  difference  between  the  desired  phase  of  quantum 
oscillations  and  the  fluctuating  phase  continuously  monitored  by  a detector  supplemented 
by  the  “calculator”  which  computes  Eqs.  (3)-(4).  The  preliminary  Monte-Carlo  results 
show  very  significant  suppression  of  decoherence  in  the  case  when  the  detector  coupling 
is  stronger  than  the  coupling  with  extra  environment. 

Several  other  predictions  related  to  experiments  with  single  quantum  systems  can  also 
be  made  using  the  Bayesian  formalism.  In  particular,  we  can  show  that  the  quantum 
oscillations  of  the  qubit  state  can  never  be  seen  directly  by  the  continuous  detector  (while 
they  can  be  computed  using  the  noisy  detector  output  and  Eqs.  (3)— (4)).  Quantitatively, 
the  spectral  peak  of  the  detector  output  I(t)  at  the  frequency  £2  = (AH1  + e2)'/2/fi  of  the 
oscillations  cannot  exceed  45/  [ ].  This  is  still  twice  as  high  as  the  classically  possible 

limit,  that  is  explained  in  the  Bayesian  approach  by  the  correlation  between  the  detector 
noise  and  the  qubit  evolution. 
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